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ABSTRACT 

We present ROSAT observations of two luminous (L 2 



10 44 ) erg s" 1 EMSS galax- 



ies, MS1019+5139 and MS1209+3917, previously classified as 'cooling flow' galaxies. 
MS1019+51 39 does not appear to be spatially extended (< f3 kpc) while its spec- 
trum is well fit by a power law with T = 1.73^q X-ray variability on a timescale 
of ~ years is also clearly detected. MSf 209+39f 7 shows no evidence of spatial exten- 
sion (< 50 kpc) but it shows variability, while its spectrum can be fit with thermal 
bremsstrahlung emission (kT = 1.8^4 kcV) or a power law model (r = 2. SOlg^ 4 ., but 
with excess photoelectric absorption above the Galactic value). All the above argue 
against thermal emission from a group of galaxies or a galaxy but in favour of an 
AGN (possibly BL Lac) interpretation. We conclude that no 'normal' galaxies with 
high X-ray luminosities have yet been detected in the EMSS survey that could be 
significant contributors to the X-ray background. 
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1 INTRODUCTION 

Although deep R OSAT surveys have shown that the major- 
ity of soft X-ray (< 2 keV) sources are type I AGN, i.e. QSOs 
and Seyfert 1 galaxies (e.g. Shanks et al. 1991, Georgan- 
topoulos et al. 1996) several galaxies have also been detected 
(Boyle et al. 1995, Carballo et al. 1995, Griffiths et al. 1996). 
Their X-ray luminosities (L x ~ ]^q42-42.5 gr g g — ij arg Qr _ 

ders of magnitude higher than those of nearby galaxies with 
similar optical magnitudes (Fabbiano 1989). Most of these 
luminous galaxies have (weak) narrow emission lines (e.g. 
Boyle et al. 1995, Griffiths et al. 1996) while few have the 
spectra of early-type galaxies, with absorption features only 
present (e.g. Griffiths et al. 1995). ft is not clear whether 
the X-ray emission is thermal or whether it emanates from 
an active nucleus. Obviously the detection of galaxies with 
such high X-ray luminosities has important implications for 
both our understanding of the physical processes in galaxies 
and the origin of the X-ray background. 

Several luminous X-ray galaxies have also been de- 
tected in the EINSTEIN Extended Medium Sensitivity Sur- 
vey (EMSS) (Gioia et al. 1990, Stocke et al. 1991). New 
optical spectra from Halpern et al. (1995) show that most 
of the objects with L x > 10 42 erg s _1 are intermediate-type 
Seyferts (e.g. Sy 1.9). Another class of luminous galaxies 
in the EMSS is that of 'cooling flow' galaxies. Stocke et 
al. classify an object as a 'cooling flow' galaxy if the op- 
tical spectrum is similar to the spectra of the dominant 
galaxies in cooling flow clusters while there is no associ- 
ated rich cluster of galaxies. Their optical spectra contain 
weak starlight absorption and strong (Wa > 25 A), low 



ionization emission with [On] 2> [Oiii]. There are at least 
four cooling flow galaxies in the EMSS, all with high X-ray 
luminosities (L x ~ iq 4 3- 44 er g s _1 ). Two of these cool- 
ing flow galaxies, MS1019+5139 and MS1209+3917, appear 
as point-like objects on Palomar E POSS plates (Maccac- 
aro et al. 1994). Their V CCD magnitudes are 18.09 and 
20.0 respectively (Stocke et al. 1991). Deep (R<23) ob- 
servations are also available (Gioia & Luppino 1994), yet 
neither source appears to be associated with a cluster of 
galaxies. These observations show that MS1019+5139 is a 
point-like source while the MS1209+3917 appears to be as- 
sociated with a faint galaxy. Initial follow-up observations 
showed that they both present [On] lines with equivalent 
widths Wa > 20A (Stocke et al. 1991). H Q was not de- 
tected in MS1019+5139 while it was out of the observed 
wavelength range in MS1209+3917. The original reported 
redshifts are z = 0.141 and z = 0.331 for MS1019+5139 and 
MS1209+3917 respectively (Stocke et al. 1991). However, af- 
ter our X-ray observations, new spectroscopic observations 
(Perlman et al. 1996a and Stocke 1996 priv. comm.) show 
that the [On] emission line in MS1019+5139 is quite weak 
(Wa = 14 A) while H a and [N n]AA6548, 6583 emission is 
found with Wa = 8 A. In the light of the new optical data, 
Perlman et al. (1996a) considered a possible BL Lac inter- 
pretation for MS1019+5139, however noting that the pres- 
ence of weak emission lines with Wa > 5A would be unique 
among X-ray selected BL Lacs. The new spectroscopic ob- 
servations (Stocke 1996 priv. comm.) fail to confirm the pres- 
ence of an oxygen emission line in MS1209+3917. Instead, 
absorption lines are detected, which give a new redshift of 
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z = 0.6 for this object; we adopt this value for the redshift of 
MS1209+3917 throughout this paper. Finally, both objects 
were detected at radio wavelengths (5 cm) with the VLA. 
Their fluxes are 2.4 and 4.6 mjy respectively (Stocke et al. 
1991). Perlman et al. (1996a) also present Position Sensitive 
Proportional Counter (PSPC) data for MS1019+5139 from 
the ROSAT all-sky survey (Molendi et al. in preparation). 
They find a hardness ratio that corresponds to a power-law 
energy index of 0.52. Unfortunately the very low exposure 
(~ 500 s) does not allow a detailed fit to the spectral data. 

Here, we discuss long exposure (~ 10 ksec) pointed 
PSPC and High Resolution Imager (HRI) observations of 
these galaxies in an attempt to constrain the origin of the 
powerful X-ray emission, and to investigate possible links 
between these and the luminous galaxies detected in deep 
ROSAT surveys. The new PSPC obsevations can provide 
useful constraints on the X-ray spectrum of these objects 
while the HRI data are important in constraining their size. 
Finally, the use of both the PSPC and the HRI data allows 
us to investigate the flux variability over long (year) and 
short (hour) timescales. 

Throughout this work, values of 50 km s _1 Mpc _1 and 
0.5 are used for Ho and qo respectively. 
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Figure 1. The HRI Radial profile of MS1019+5139 (with 1 a 
errors). The profile of AR Lac is shown (solid line). A 99 per cent 
upper limit to the linear size of the object was found to be 12.8 kpc 
(see text for details). Contributions from each bin to the x 2 (from 
a comparison of the AR Lac data and the MS1019+5139 data) 
are shown in the lower panel. 



2 THE ROSAT OBSERVATIONS 

MS1019+5139 was observed with the PSPC detector (Pf- 
effermann et al. 1986) on board ROSAT (Trimmer 1983) 
on 21st October 1993 while MS1209+3917 was observed 
on the 3rd June 1994. After rejecting the data for times 
of high particle rates (Master Veto rate > 170, e.g. Snow- 
den et al. 1992) we end up with exposure times of 7488 
and 11857 s for MS1019+5139 and MS1209+3917 respec- 
tively. For the estimated fluxes and luminosities, see tables [j] 
and |2[ Note that the fluxes and luminosities quoted through- 
out the paper are not corrected for absorption. The source 
centroids were estimated using the intensity weighted cen- 
troid algorithm in the Starlink Asterix package. The es- 
timated celestial coordinates (J2000) are a = 10 h 22 m 12!4, 
5 = +51°24'1'.'0 for MS1019+5139 and a = 12 h ll m 34:3, 
5 = +39°00 / 55'.'0 for MS1209+3917, with an error circle 
radius of ~ 5 arcsecs. The proposed optical counterparts 
are offset from these positions, in the R.A. and dec. direc- 
tions, by +2.3, -0.7 and -1.3, -0.7 arcsec for MS1019+5139 
and MS1209+3917 respectively. Positive offsets indicate op- 
tical positions to the North and East of the X-ray centroids. 
A PSPC serendipitous observation of MS1209+3917 is also 
used (ROR rp 700725), obtained on 31st May 1991, having 
an exposure time of 38125 s, with MS1209+3917 at a po- 
sition ~ 22 arcmin off-axis but not obscured by the PSPC 
window support structure. These data were retrieved from 
the University of Leicester LEDAS ROSAT public archive. 

Data from two ROSAT HRI pointed observations are 
used, firstly an on-axis observation of MS1019+5139 (ROR 
rh600125) with an exposure of 10924 s, taken between 26th 
October 1991 and 18th July 1992 and an observation of 
MS1209+3917 obtained on 18th May 1995 with an expo- 
sure of 11360 s (ROR rh701844). 



3 RESULTS 

3.1 MS1019+5139 

In Fig. 1 we present the radial count distribution for 
MS1019+5139 from the HRI observation. Data were binned 
into 1 arcsec bins, out to, and including, a radius of 18 arc- 
sec, the radius which encloses w 85 per cent of the source 
photons. To determine an upper limit to the angular size of 
the object, we modelled the data as a simple Gaussian dis- 
tribution of variable sigma. An upper limit to the gaussian a 
is found, at the 99 per cent confidence level, to be 3.45" , cor- 
responding to a linear distance, at the redshift of the object 
(z = 0.141), of 12.8 kpc. For comparison, the radial count 
distribution of an on-axis HRI observation of AR Lac (a 
point source used in calibrating the ROSAT HRI) is plotted 
as a solid line. A x 2 value was found from comparison of the 
MS1019+5139 HRI radial profile and that of AR Lac (nor- 
malizing the background subtracted data from AR Lac to 
the background subtracted data from MS1019+5139), giv- 
ing xt = 1-06 from 17 degrees of freedom. 

We have divided the PSPC and the HRI data into in- 
tervals of 1 hour and 400 s respectively, in order to look 
for short-term variability. The results are shown in Fig. 2. 
It can be seen there is no evidence for short-term variabil- 
ity, however, comparison of the PSPC flux (2.75 ±0.10 x 
10" 12 erg cm- 2 s _1 ) with the HRI fluxes (2.76 ± 0.26 x 
10~ 12 erg cm" 2 s _1 and 4.88±0.32 x 10~ 12 erg cm"V) and 
the ROSAT all-sky survey flux from Perlman et al. (1996a) 
(4.03 x 10~ 12 erg cm~ 2 s~ 1 ) suggests variability on the long- 
term (year) timescale. The above fluxes refer to the 0.1-2.0 
keV band, and were estimated using the best fit power-law 
model described below. A \ 2 t es t rejects the null hypothesis 
of no long-term variability at the 99.99 per cent confidence 
level. We note that the observed variability cannot be due 
to uncertainties or variations in the calibration of the PSPC 
and the HRI. Hasinger et al. (1993) present the results of 
a long term program of calibration observations with the 
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Figure 2. The variation of source flux with time for 
MS1019+5139, from both the PSPC observation (right hand 
panel) and the two parts of the HRI observation (left and cen- 
ter panels). The mean fluxes are shown as dotted lines and the 
PSPC all-sky survey flux from Perlman et al. (1996a) is denoted 
by the arrow (see text for details). Fluxes were found in the en- 
ergy range 0.1-2.0 keV, and used the best fit power-law model 
(r = 1.73, N H = 1.34 X 10 2() cm- 2 ). The x-axis shows the time 
from the start of each separate observation, in hours. Errors on 
the data are 90 per cent errors, from photon statistics and uncer- 
tainties in the model used for converting count rates to fluxes. 



PSPC, concluding that over a 3 year period, the count rate 
of the supernova remnant N132D measured in the PSPC 
was constant within 2 per cent. Similar work using the HRI 
(David et al. 1996) finds that the observed count rate of 
N132D to vary by ±5 per cent over the period 1990-1995. 
We also note that there is an excellent agreement (better 
than 1 per cent) between the HRI and PSPC flux calibra- 
tion. This is derived by comparing the PSPC and HRI flux 
of the galaxy cluster A2256 (David et al. 1996), where the 
count rates are converted to flux using the X-ray spectrum 
from Markevitch & Vikhlinin (1997). In contrast, the count 
rate of MS1019+5139 in the HRI increases by 70 per cent in 
~ 6 months and a comparison of the two PSPC fluxes shows 
a 40 per cent decrease between the all-sky survey and the 
pointed observation. We choose not to include in our analy- 
sis the Einstein IPC measurement obtained about 10 years 
earlier, due to differences in the calibration of the responses 
of the two satellites (Fiore et al. 1994). 

Spectral fits were made with the simplest possible mod- 
els: a power law and a Raymond-Smith hot plasma model 
(Raymond & Smith 1977), both with an absorption com- 
ponent included, having the column density either fixed to 
the Galactic value, or as a free parameter. The results are 
presented in table 1. We see that the power law model 
provides a marginally better fit compared to the plasma 
model, albeit not statistically significant. The value for the 
F-statistic is 0.818, implying that a power-law model gives 
a better fit only at the ~ 30 per cent level. The derived 
value for the hydrogen column is Nh = 1.34 x 10 20 cm~ 2 , 
with r = 1.73, somewhat above, but within the 90 per cent 
errors of the Galactic Nh value from Stark et al. (1992) of 
Nh = 0.88 x 10 20 cm~ 2 . The photon index is consistent, 
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Figure 3. The spectral data for MS1019+5139, shown with the 
best fitting power-law model (solid line) and the residuals from 
this model plotted below. 



within ~ 2cr errors, with the photon index from Perlman et 
al. (1996a), of F = 1.52. The folded spectrum for the best fit 
power law model, along with the residuals, is given in Fig. 
3. We find that the best fit temperature of the Raymond- 
Smith model is kT = 3.7 keV using the Galactic absorption 
column density and the heavy element abundance fixed at 
0.3 solar units. Allowing the heavy element abundance to be 
a free parameter, we find a best-fit value of 0.0 solar units, 
the 90 per cent upper limit being 0.2 solar units. The derived 
temperature is not significantly changed. 



3.2 MS1209+3917 

The observed radial count distribution from the HRI data is 
given in Fig. 4, together with the radial count distribution 
found from an on-axis observation of AR Lac. Comparison 
of the MS1209+3917 data with the AR Lac data gives a 
xt value of 0.8 for 17 degrees of freedom, clearly showing 
that our object is a point source. Using a similar method 
to determine the radial extension as for MS1019+5139, we 
found an upper limit of a = 4.06", giving a value for the 99 
per cent upper limit to the linear size of the object of 49.5 
kpc. 

Like MS1019+5139, this source is seen to vary on 
year timescales, (Fig. 5), evidenced by the fluxes found, 
in the 0.1-2.0 keV band, from the pointed PSPC obser- 
vation, the serendipitous observation and the HRI data, 

\8.69 ± 1.61 x 



which are 6.73 ± 0.35 x 10~ 13 erg cm~ 2 



10~ 13 erg crrrV 1 



and 3.43 ±0.31 x 10" 13 erg cm" 2 s _1 re- 



spectively. The conversion from count-rate to flux was per- 
formed using the best-fit thermal model, described below. 
No significant variation in flux is seen throughout the dura- 
tion of each individual observation. 

Spectral fitting was performed using first a power law 
model and then a Raymond-Smith plasma model. The re- 
sults are presented in table 2, and Fig. 6 shows the folded 
spectrum with the best fitting plasma model. 

Both models give acceptable fits when the column den- 
sity is a free parameter, although the plasma model does not 
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Figure 4. The HRI radial profile of MS1209+3917. Data are 
shown with 1 a errors. The profile of AR Lac is shown as a solid 
line. A 99 per cent upper limit to the linear size of the object 
was found to be 49.5 kpc (see text for details). In the lower panel 
are the contributions to the \ 2 from each bin, caculatcd from a 
comparison of the AR Lac data with the MS1209+3917 data. 
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Figure 5. The variation with time of the flux of MS1209+3917, 
from the serendipitous PSPC observation (left hand panel), the 
pointed PSPC observation (center panel) and the pointed HRI 
observation (right hand panel). The mean fluxes are shown as 
dotted lines. Count rates were converted to fluxes in the energy 
range 0.1-2.0 keV using the best fit Raymond-Smith model. 



statistically require an absorption column greater than the 
galactic value of N H = 1.86 x 10 2() cm" 2 , (Stark et al. 1992). 
However, we find the column density derived for the power- 
law model is substantially higher (Nh = 4.85 x 10 20 cm -2 ) 
than the observed Galactic value, and this excess above the 
Stark et al. value is significant at over the 90 per cent level. 
The corresponding photon index is T = 2.50. When we fix 
the Nh value to the Stark et al. value we obtain a value of 
r = 1.74 for the power-law index. The best fit temperature 
found with the plasma model with fixed column density and 
the heavy element abundance fixed at 0.3 solar units is 1.8 
keV. The heavy element abundance is not well constrained: 



the best-fit value is 0.56 solar units with an upper and lower 
limit of 1.8 and 0.2 solar respectively at the 90 per cent 
confidence level. 



4 DISCUSSION 
4.1 MS1019+5139 

The detection of variability in the R OS AT observations im- 
mediately discriminates against a thermal origin for the ra- 
diation, ruling out emission from hot diffuse gas in an early- 
type galaxy or a group of galaxies. A hot intergalactic gas 
origin is further ruled out by the lack of spatial extension in 
the HRI observations which limits the size of the source to 
~ 13 kpc. Hence, the X-ray data favour an AGN interpre- 
tation. 

The absence of strong emission lines in the optical spec- 
trum (weak [O n] and H a + [N n] only are present) argues 
against a QSO or Seyfert source and could suggest possi- 
bly a galaxy similar to those found in deep R OSAT surveys, 
or alternatively a BL Lac object. We note however that, 
the X-ray luminosity of MS1019+5139 (L x > 10 44 erg s" 1 ) 
is around an order of magnitude brighter than the galax- 
ies found in deep ROSAT surveys. Perlman et al. (1996a) 
first suggested a BL Lac identification based mainly on the 
optical spectra and the lack of spatial extension in X-rays. 
However, they noted that the equivalent widths of the emis- 
sion lines are larger than those of X-ray selected BL Lacs. 
We further note that the X-ray spectrum derived here is 
well represented by a power-law with F = 1.73lg']J, much 
flatter than the canonical X-ray selected BL Lac spectrum 
of T ~ 2.52 (Padovani & Giommi 1996). The detection of 
MS1019+5139 in radio wavelengths adds support to a BL 
Lac scenario, as most X-ray selected BL Lacs show strong 
radio emission (Perlman et al. 1996b). Indeed, with values 
of 1.06 and 0.25 for a ox and a ro respectively (Stocke et al. 
1991), the location of this object on the a ox — «ro plane (op- 
tical to X-ray vs. radio to optical spectral index) is very close 
to the X-ray selected BL Lac region, approximately delim- 
ited by the values 0.6 < a ox < 1.2 and 0.3 < a ro < 0.5 (see 
Fig. 6 of Stocke et al. 1991). We emphasize however that the 
X-ray data alone cannot provide a definitive identification 
for this object. Further polarization and radio observations 
are necessary to prove that MS1019+5139 is a BL Lac ob- 
ject. 



4.2 MS1209+3917 

The lack of X-ray extension strongly argues against ther- 
mal emission from intergalactic gas in a group or cluster 
of galaxies. The size of MS1209+3917 is <50 kpc at the 
99 per cent confidence level, compared to ~260 kpc for 
the cluster 0055-279 (Roche et al. 1995), which is at sim- 
ilar distance (z=0.561) and has a comparable luminosity, 
L x « 10 44 erg s _1 . Furthermore, MS1209+3917 does not 
show any evidence for a cluster on deep CCD images in clear 
contrast with the 0055-279 case. Thermal emission from a 
compact group of galaxies is unlikely as the X-ray lumi- 
nosity of MS1209+3917 is that of rich clusters of galaxies, 
while typical group luminosities are L x < 10 43 erg s -1 (e.g. 
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Figure 6. The observed spectral data for MS1209+3917, with 
the best fit Raymond-Smith model and residuals. 

Ebeling, Voges, Bohringer 1994). Moreover, the strongest ar- 
gument against thermal emission comes from the detection 
of long-term (year timescale) variability between the PSPC 
and the HRI observation. 

The X-ray luminosity of MS1209+3917 is approxi- 
mately two orders of magnitude higher than that of the 
most luminous nearby, optical galaxies observed by EIN- 
STEIN (Fabbiano 1989) and about an order of magnitude 
higher than the galaxies detected in deep ROSAT surveys 
(Griffiths et al. 1995, Moran et al. 1996). All the above 
strongly suggest the presence of an active nucleus. The 
optical spectrum, which is unfortunately unpublished yet, 
holds important clues on the nature of this object. The 
absence of emission lines makes a BL Lac intepretation 
a strong possibility. The X-ray spectrum power-law pho- 
ton index (r = 2.501q;42) i s similar to the average 0.1-2.0 
keV photon index of X-ray selected BL Lacs (r = 2.52) 
(Padovani & Giommi 1996). We observe significant absorp- 
tion (N H = 4.85ti' ( 62 x 10 20 cm" 2 ) above the Galactic value 
of Nh = 1.86 x 10 20 cm~ 2 ; such an excess absorption is not 
atypical in BL Lacs (Urry et al. 1996). The BL Lac hypoth- 
esis is further supported by the broad band spectral indices 
(q ox = 1.05, Q ro = 0.44), which place MS1209+3917 well 
within the X-ray selected BL Lac region on the a ox — a TO 
plane. 
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Table 1. Spectral fitting results for MS1019+5139. The ab- 
sorption columns denoted by a dagger are the Galactic values 
from Stark et al. (1992) and arc in units of (10 20 cm -2 ). F x is 
the flux in the energy band 0.1-2.0 keV, expressed in units of 
10~ 12 erg cm~ 2 s _1 . L x is the luminosity, in the same energy 
band, expressed in units of 10 44 erg s~ 1 . The plasma tempera- 
ture, kT, is in keV. 
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Table 2. Spectral fitting results for MS1209+3917. The ab- 
sorption columns denoted by a dagger are the Galactic values 
from Stark et al. (1992) and are in units of (10 2() cm -2 ). F x is 
the flux in the energy band 0.1-2.0 keV, expressed in units of 
10 — 13 erg cm _2 s _1 . L x is the luminosity, in the same energy 
band, expressed in units of 10 44 erg s — 1 . The plasma tempera- 
ture, kT, is in keV. 



from a 'normal' galaxy or a cluster of galaxies and clearly 
support an AGN interpretation. The possibility that there 
is some connection to the non-QSO faint sources detected 
in deep ROSAT surveys (e.g. Boyle et al. 1995) cannot be 
definitively ruled out on the basis of the X-ray data alone. 
However, the absence of strong, broad emission lines in the 
optical spectra as well as the detection of radio emission, 
possibly favours a BL Lac interpretation. Of course, further 
variability and polarization observations in the optical and 
radio wavebands are necessary to confirm the BL Lac hy- 
pothesis. We note that if MS1019+5139 is found indeed to 
be a BL Lac then it will be an atypical one both because 
of its flat X-ray spectrum and its optical emission lines and 
will be worthy of further study. 



5 CONCLUSIONS 

We have presented the results of ROSAT observations of 
two X-ray highly luminous galaxies (MS1019+5139 and 
MS1209+3917) from the EMSS, previously classified as 
'cooling flow' galaxies, in an attempt to constrain the ori- 
gin of the powerful X-ray emission. The HRI observations 
suggest that neither of the two objects shows evidence of 
spatial extension. The PSPC spectra are consistent with a 
power-law model, although a bremsstrahlung model cannot 
be strongly ruled out. Finally, long-time variability is clearly 
detected between the PSPC and HRI observations in both 
objects. The above together with the high X-ray luminosities 
of these objects strongly rule out a thermal emission model 



In summary, the ROSAT observations provide pow- 
erful constraints on the nature of these objects. Although 
the exact identification remains evasive, we are able to rule 
out thermal emission models from a 'normal' galaxy i.e. a 
galaxy without an active nucleus. We conclude that no 'nor- 
mal' galaxy has been detected at high X-ray luminosities, 
L x > 10 43 erg s~\ in the EMSS or elsewhere, extending 
previous results by Halpern et al. (1995). 
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